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Abstract 

The COMPASS collaboration at CERN has measured diffractive dissociation of 190GeV/c 
pions into the 7r“7r“7r"'' final state using a stationary hydrogen target. A partial-wave analysis 
(PWA) was performed in bins of Stt mass and four-momentum transfer using the isobar model 
and the so far largest PWA model consisting of 88 waves. A narrow = 1++ signal is 
observed in the /o(980) tt channel. We present a resonance-model study of a subset of the 
spin-density matrix selecting Stt states with = 2++ and 4++ decaying into p{770) tt 
and with = 1++ decaying into /o(980)7r. We identify a new oi meson with mass 
(1414^i[3)MeV/c^ and width (153^23)M®V/c^. Within the final states investigated in our 
analysis, we observe the new ai(1420) decaying only into /o(980)7r, suggesting its exotic 
nature. 
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One of the great challenges in particle physics is the understanding of how hadronic matter is 
constructed from its basic building blocks, quarks and gluons. Although many possibilities are 
allowed, which all follow the principle of confinement, almost all hadrons observed can be described 
by the constituent quark model. The known light-meson spectrum is presently interpreted in 
terms of qq quark-model states that are associated with flavor SU(3) multiplets according to their 
mass and quantum numbers. For some combinations, more states were reported than 
can be accommodated by SU(3) symmetry. Depending on their coupling to specific production 
mechanisms and their decay pattern, these states are interpreted as either carrying a strong 
glueball component, e.g. /o(1500), as molecular-type excitations, e.g. /i(1420), or as tetra-quark 
states. For a detailed review see Ref. [T]. However, their exotic structure is under debate, unlike 
for states that carry spin-exotic quantum numbers, e.g. = 1 and hence cannot be qq 
states. This is in contrast to the sector of heavy mesons containing c or 6 quarks, where exotic 
mesons have clearly been identified, e.g. X, Y, and Z-states. In particular, the recent observation 
of charged Z-states, such as Zc(3900)^ and Zf,(10 610)^, has proven the existence of mesons 
with exotic structure EElIl]. Already the existence of one system with a wave function and/or 
quantum numbers requiring an explanation beyond “ordinary” mesons implies the existence 
of a large number of additional matter states. Unless new principles for the construction of 
color-neutral hadrons are found, this should hold for all quark-flavour combinations. 

In the sector of light mesons, the issue of exotic states remains unresolved. The lowest-mass 
state discussed in this context is the /o(980), which contains nn (n = {u,d}) and a dominant 
ss component. It has also been interpreted as a KK molecule 0 El E]- The /i(1420) with a 
width of only 55MeV j(? couples strongly to KK and was also suggested to be a molecular-type 
structure j8]. In Ref. [3, the Particle Data Group has tentatively adopted the scenario of /i(1420) 
being the SU(3) partner of /i(1285). In the class of spin-exotic mesons, the 7ri(1600) is the only 
meson observed by several experiments in different decay modes. However, the masses quoted 
and in particular the widths vary considerably between different experiments, and values for the 
width often exceed 400 MeV/. This suggests the existence of dynamical effects similar to the 
case of ai(1260). The situation is characterized by individual states without recognizable pattern 
and, except for ao(980) and /o(980), the absence of any observed isospin partners. 

In order to search for new exotic mesons, we have analyzed the diffractive reaction tt~ -|- p —)• 
7r“7r“7r+ + Precoil with the focus on waves with quark-model combinations]^ We have studied 
the states in order to search for a possible partner of the isosinglet /i(1420). Our 

analysis aims at the charged isospin 1=1 analogue decaying into 7r“7r“7r'’'. Although this final 
state and the mass range of 1 to 2GeV j(? have already been studied by many experiments, see 
e.g. Refs. |10l ITTl I12| . the improvement by almost two orders of magnitude in sample size has 
opened a new avenue for analysis. 

The COMPASS experiment |T3l [T3] is located at the M2 beam line of the CERN Super Proton 
Synchrotron. For this measurement, we used a beam of 190GeV/c 7r“ with 96.8% purity, 
impinging on a 40 cm long liquid-hydrogen target that was surrounded by a recoil-proton detector 
(RPD). Incoming pions were identihed with a Cherenkov counter placed in the beam line at 
the entrance to the experimental area. The large-acceptance high-precision spectrometer is well 
suited for investigating high-energy reactions at low to intermediate values of which denotes 
the reduced squared four-momentum transfer to the target proton with t' = |t| — |t|min- For this 
measurement, t' is chosen to be in the range from 0.1 to 1.0 (GeV/c)^, where the lower bound is 
dictated by the acceptance of the RPD and the upper bound by the steep drop of the number of 
events with increasing t'. Outgoing charged particles are detected by the tracking system and 
their momenta are determined using two large-aperture magnets. 

^The C-parity refers to the neutral state of the isospin multiplet. 
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The data presented in this Letter were recorded in the year 2008. A detailed description of setnp, 
data selection, and analysis scheme can be found in Refs. |15[ I16| . The trigger is based on a 
recoil-proton signal in the RPD in coincidence with an incoming beam particle and no signal 
in the beam-veto counters. We require a production vertex located within the target volume, 
with one incoming beam-pion track and three outgoing charged particles, compatible with the 
pion hypothesis based on information from the RICH counter. The sum of the momenta of the 
outgoing particles is required to be equal to the average beam momentum within two standard 
deviations, i.e. ±3.78GeV/c. We require Feynman-x of the fastest final-state pion to be below 0.9 
for rapidity differences between the fast 7r“ and the slower 7r“7r+ pair in the range 2.7 to 4.5. 
This suppresses the small contamination of centrally produced final states, which contribute 
mainly at higher Stt masses. A total of 46 x 10® events was selected in the mass range between 
0.5 and 2.5 GeV /(?. 


In order to extract the spectrum of resonances produced in the reaction, we have performed a 
partial-wave analysis (PWA) that is pursued in two steps. First, we fit the intensity distributions 
in the 5-dimensional phase space independently in one hundred 20 MeV/wide bins of Svr mass 
msTT, each divided into 11 bins of t'. We adopt the notation M'^ [isobar] it L to define partial 
waves. Here, £ denotes the reflectivity and M > 0 the magnitude of the spin projection along the 
beam axis (see Ref. CZl), while L is the orbital angular momentum between the isobar and the 
bachelor pion in the decay of the Svr state. We use the isobar model, which for our hts contains 
88 waves, i.e. 80 waves with positive and 7 with negative reflectivity as well as one non-interfering 
wave representing three uncorrelated pions. This set contains all significant isobars that decay 
into 7r“7r^ and has been derived from a larger set with 128 waves by requiring a minimum relative 
intensity of about 10“^. The likelihood ht function is built from two incoherently added terms 
that correspond to the two values of reflectivity, e = ±1. Each term coherently sums over all 
partial-wave amplitudes that belong to the respective value of e. Details on the fit model, the 
fitting procedure, and the results are described in Refs. |15| I16j . The division of our data set 
into 11 bins of t' is motivated by the very different t'-dependences of resonant and non-resonant 
components [lOl I15| . In all partial waves studied, the intensity of non-resonant, i.e. Deck-like 
components |18) . drops off much faster with increasing than that of resonances. 


The result of this first PWA step reveals a number of well-known resonances with = 0 2 

2+'*', and They are identihed by structures in the mass spectra and a mass-dependent 

phase that is measured against the reference wave T*”*" O"*"/9(770) vr S'. The 1 0^/o (980) tt P 

intensity shows a clear signal slightly above 1.4 GeV/c^, which cannot be associated with a known 
1state [see data points in Fig. |l(a)| . Rapid phase rotations with respect to known resonances 
are observed in the signal region, independent of t' [see points in Figs. |l(d)| and l(e)| . The same 
feature is observed in the hnal state 1191. 


In the second analysis step, we use a resonance model to ht the spin-density matrices resulting 
from the hrst analysis step simultaneously in all bins of t' and in wave-specihc ranges in m 3 ,r- 
Typically, only subsets of these spin-density matrices are ht simultaneously. In this Letter, we 
present such a ht using a minimal set of 3 waves, namely 1"*“ p{770) vr D, 1^ p(770) vr G, 
and 1++ O'*' /o(980) n P. The hrst two waves contain the known 02(1320) and 04(2040). These 
two waves act as interferometers in order to search for structures in I"*"*" 0^ /o(980) ttP, where no 
resonances have yet been reported. For this ht, we model the amplitude of each wave by a coherent 
superposition of a resonant contribution, which is described by a relativistic Breit-Wigner (BW) 
amplitude and a non-resonant contribution. Hence both components are allowed to interfere. In 
the 4'^'^ and 1+'’' waves, the latter are described by terms of the form P'(m 3 ,r) = 
where ci is a ht parameter and q is the two-body break-up momentum for a particular isobar 
at the mass For the non-resonant term in the 2'^~^ wave, this parametrization is extended 
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Figure 1: (Color online) Results of the PWA in Stt mass bins of 20MeV/c2 width (data points with 
statistical errors only). The red curves in the intensity distributions [panels (a)-(c)] represent the fit 
model, which is the coherent sum of resonances (blue) and non-resonant contributions (green). The fit 
is constrained to the mass range indicated by the continuous curves. Extrapolations of the model are 
shown as dashed curves. Panel (d) shows the relative phase between 1++ and 4++ together with the 
model curves and (e) the phase between two 1++ decay modes. Here, the different colors distinguish 
3 exemplary t' bins. The t'-dependence of the ai(1420) intensity is shown in (f) together with the result 
of a single-exponential fit (red line) yielding a slope parameter of 6 « 10 (GeV/c)^^. 


to include an explicit t'-dependence. A simple BW amplitude is used for the 04(2040) and the 
= I"'"'' state, a BW amplitude with mass-dependent width for the 02(1320). In the latter, 
the decay phase space is approximated assuming quasi-two-body decays into 80% p{770) vr and 
20 % 77 vr. 


The result of this fit is shown as curves in Figs. 1(a) to |l(d)[ in which the model curves describe 
the data well. The blue curve in Fig. 1(a) reveals the existence of a new axial-vector state in 
the I"*”*" 0^ /o(980) tt P wave, which we introduce as oi(1420). This wave collects only 0.25% of 
the total observed intensity. Its resonance interpretation is supported by the observation of a 
rapid phase variation by about 180° across the peak region with respect to the 4'^^ [see Fig. l(d)| 
and 2~^~^ reference waves. As illustrated in Fig. 1(e), the I"*"*' 0"*“ /o(980) tt P wave shows similarly 
rapid phase motion also relative to the 1"’“"’' 0"^ p(770) vr S wave. This indicates that the observed 
structure in the /o(980) vr decay mode is not caused by the high-mass tails of the ai(1260), which 
dominates the p{770) tt wave. Our fit reveals a BW mass for the ai(1420) of 1414 MeV j(? and 
a width of 153 MeV j(?. The observed shift of the measured peak position with respect to the 
resonance position in Fig. 1(a) is due to destructive interference of the BW resonance (blue) and 
the non-resonant term (green). In this wave, the fit model is chosen to cover the mass range up to 
1.6GeV/c^. The observed tail of the intensity at higher masses may be attributed to non-resonant 
contributions that are not described by the present model. 
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The resonance-model fit is performed simultaneously in all 11 bins of t'. We allow production 
strengths and phases of resonances and non-resonant contributions to vary with t'. Spectral 
shapes and BW parameters are assumed to be independent of t'. Figure [l (f) | shows the resulting 
t'-spectrum of the production intensity of the BW amplitude that describes the ai(1420). The 
BW intensity and that of the non-resonant contribution show a steep, approximately exponential 
f'-dependence. Fits with a single exponential yield the slope parameters. Resonances are typically 
described by slope parameters 6 ~ 8 to 10(GeV/c)^^ that are smaller than those of the non¬ 
resonant contributions with 6 ~ 12 to 15 (GeV/c)^^ [15| . The new ai(1420) has a slope parameter 
of 6 ~ 10 (GeV/c)^^ that is similar to those of 02(1320) and 04(2040), which supports its resonance 
interpretation. The fact that the Oi(1420) is produced with nearly constant phase offset relative 
to the 02 ( 1320 ), independent of t', provides further support for this interpretation. As expected, 
the slope of the non-resonant contribution in the l'’'+ wave is steeper with 6 ~ 13 (GeV/c)^^. 

The 88 partial-wave set contains three independent contributions for the tttt S'-wave isobars, 
namely the /o(980) with parameters taken from Ref. |20) . a broad component denoted [7r7r]5 taken 
from elastic vrvr 5-wave scattering |21) . and the /o(1500) described by a simple BW amplitude. 
The Oi(1420) is observed only in l^^ 0^ /o(980) ttP, while no signal with corresponding phase 
motion is seen in I"*”*" O"*" [vrvrjs'vr P or any other I"*”*" wave. In order to confirm the unique 
coupling of Oi(1420) to /o(980) tt, we have investigated in a separate study [15[ [22] the structure 
of the 7r“7r"’“ subsystem forming 0"’“'*' isobars using a novel fit procedure for the partial-wave 
decomposition in bins of and f. Instead of describing the O"*"*' isobars by several amplitudes 
with fixed shape, their mass dependence is replaced by a piecewise constant function across 7712 , 7 , 
which is determined from data. We thus remove possible bias originating from the isobar model 
for O"’"*'. For 7773 ,, around the new resonance, a clear intensity correlation of the new ai(1420) 
with the /o(980) is observed within the extracted 0"’“"’' isobar amplitude |15) . 

Due to the large data set, statistical uncertainties of the extracted resonance parameters are 
negligible compared to systematic ones. Therefore, we performed extensive systematic studies 
concerning event-selection cuts, the model used in the first step of the PWA fit, as well as 
wave set and parametrizations employed in the resonance-model fit. The result is stable under 
all these studies. The main systematic uncertainties arise from the resonance-model fit. The 
estimated total systematic uncertainty is ^jgMeV/c^ for the Oi(1420) mass and for 

the width. This was obtained by changing the set of waves used in the resonance-model ht, the 
event selection criteria, and the parametrization of the non-resonant terms. Using instead of 
a simple BW amplitude one with a mass-dependent width, with /o(980)7r as the only decay 
channel, yields central values for mass and width of 1433 MeV j(? and 146 MeV/c^, respectively. 
Pseudo data simulated using the full 88 waves with and without the l"*"^ 0"*“ /o(980) tt P wave, 
which were subsequently analyzed with our standard 88-wave PWA, show no indication for model 
leakage artihcially populating the l"*"^ O"*" /o(980) vrP wave. 

We have estimated the statistical signihcance of the ai(1420) signal by rescaling the statistical 
uncertainties of the data points such that the fit has a y^/ndf = 1 (following Ref. [9]). The 
decrease of the probability of a ht performed without the assumption of a new resonance in 
the I"*”*" wave translates into a signihcance of the ai(1420) measurement, which by far exceeds 
the value of 5 a. The same result is obtained when only those spin-density matrix elements are 
included in the calculation, which contain the l'’'+ wave. Therefore, the result does not depend 
on possible imperfections in the description of the 2~^~^ and reference waves at masses much 
below or above the new ai(1420). 

Summarizing our analysis, we have performed a resonance-model ht based on a spin-density 
submatrix that was obtained by the so far most extensive 88-wave Svr PWA using the large GOM- 
PASS data set. Restricting this resonance-model ht to the three waves I"*" p{770) tt D, I"*" 
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p(770) vr G, and O"*" /o(980) vr P, we have observed a new oi meson at m = (1414MeV/ 
with a width of -T = (153^23) MeV/c^. This finding was made possible by onr large event sample 
and, when compared to previons experiments, also by the more homogeneons acceptance of the 
COMPASS setnp in the five phase-space variables nsed in the PWA nans]- 

The interpretation of this new state is still unclear. Quark-model calculations including tetra- 
quark states predict nnrm and nsns iso-vector states at 1381 and 1530MeV/c^, respectively |23| . 
The molecular model for the /i(1420) proposed in Ref. jSj could possibly be extended to the 
isovector case. After our first announcement of the ai(1420) |24) . several explanations were 
proposed |25l EH EH EH], none of them describing all features of the data. The properties of the 
ai(1420) suggest it to be the isospin partner of the /i(1420). This is supported by its mass value 
of 1414 MeV !(? and by its strong coupling to /o(980), which is interpretable as a KK molecule. 
The ai(1420) width of only 153MeV/c^ is narrow as compared to most other known iso-vector 
states, which typically have widths between 250 and 350 MeV jc?. The much smaller width of the 
/i(1420) of only (54.9 ± 2.6) MeV/c^ can be explained by its strong coupling to KK with the 
corresponding phase space being much smaller than that for ai(1420) decaying into /o(980) vr. 
The Oi(1420) and the /i(1420) may possibly be the first observed isospin partners for a ttKK 
molecular-type excitation, which obey isospin symmetry. This interpretation suggests further 
experimental and theoretical studies of the ttKK final state. 
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